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We perform a statistical weak lensing analysis of dark matter profiles around tracers of halo mass 
from galactic- to cluster-size halos. In this analysis we use 170 640 isolated ~ L» galaxies split 
into ellipticals and spirals, 38 236 groups traced by isolated spectroscopic Luminous Red Galaxies 
(LRGs) and 13 823 MaxBCG clusters from the Sloan Digital Sky Survey (SDSS) covering a wide 
range of richness. Together these three samples allow a determination of the density profiles of dark 
matter halos over three orders of magnitude in mass, from 1O 12 M0 to 1O 15 M0. The resulting lensing 
signal is consistent with an NFW or Einasto profile on scales outside the central region. In the inner 
regions, uncertainty in modeling of the proper identification of the halo center and inclusion of 
baryonic effects from the central galaxy make the comparison less reliable. We find that the NFW 
concentration parameter C2006 decreases with halo mass, from around 10 for galactic halos to 4 for 
cluster halos. Assuming its dependence on halo mass in the form of C2006 = Co (M/lO 14 /i- 1 M )-' 3 
we find Co = 4.6 ± 0.7 (at z = 0.22) and /3 = 0.13 ± 0.07, with very similar results for the Einasto 
profile. The slope (/?) is in agreement with theoretical predictions, while the amplitude is about 
two standard deviations below the predictions for this mass and redshift, but we note that the 
published values in the literature differ at a level of 10-20% and that for a proper comparison our 
analysis should be repeated in simulations. We compare our results to other recent determinations, 
some of which find significantly higher concentrations. We discuss the implications of our results 
for the baryonic effects on the shear power spectrum: since these are expected to increase the halo 
concentration, the fact that we see no evidence of high concentrations on scales above 20% of the 
virial radius suggests that baryonic effects are limited to small scales, and are not a significant source 
of uncertainty for the current weak lensing measurements of the dark matter power spectrum. 



PACS numbers: 

I. INTRODUCTION 



The density profile of dark matter (DM) halos is one 
of the fundamental predictions of cosmological models 
in the nonlinear regime, determined using N-body sim- 
ulations 0, a a is, a 0, a a m hi m, m m m . 

The profile is often parametrized with the so called NFW 
profile T2|, a broken power-law characterized by a scale 
radius where the slope is approximately —2, and by the 
virial radius which parametrizes its mass. The ratio of 
the latter to the former is called the concentration c, and 
is a measure of the density of the halo in the inner regions: 
a higher concentration implies a higher density of the halo 
at a fixed fraction of virial radius. The concentration 
is predicted to be mildly dependent on mass such that 
higher mass halos are less concentrated than lower mass 
halos. Because the normalization of the concentration- 
mass relation depends on the matter power spectrum nor- 
malization, its shape, and the matter density fl m , mea- 
suring the halo profile as a function of halo mass can 
teach us about the underlying cosmology. 
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Many different methods have been used to measure the 
halo profile of clusters. Kinematic tracers such as satellite 
galaxies, in combination with a Jeans analysis or caustics 
analysis, can give information over a wide range of physi- 
cal scales and halo masses. While the issues of relaxation, 
velocity bias, anisotropy of the orbits and interlopers con- 
tinue to be debated and need to be carefully addressed, 
recent results suggest a good agreement with theoretical 
predictions [lj| Qj], [H, H3, EH . Hydrostatic analyses 
of X-ray intensity profiles of clusters use X-ray intensity 
and temperature as a function of radius to reconstruct 
the density profile. They have the benefit of thermal gas 
pressure being isotropic, but may be biased due to the 
possible presence of other sources of pressure support, 
such as turbulence, cosmic rays or magnetic fields. These 
cannot be strongly constrained for typical clusters with 
present X-ray data [22j , but could modify the hydrostatic 
equilibrium and affect the conclusions of such analyses. 
Recent results are encouraging and are in a broad agree- 
ment with predictions, although most require concentra- 
tions that are hi ghe r than those predicted by a concor- 
dance cosmology [23l. l24l. [25| . While the above-mentioned 
systematic biases cannot be excluded, the small discrep- 
ancy could also be due to baryonic effects in the central 
regions, due to selection of relaxed clusters that may be 
more concentrated than average [2g, or due to the fact 
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that at a given X-ray flux limit, the more concentrated 
clusters near the limi ting mass are more likely to be in- 
cluded in the sample [2q . 

Gravitational lensing is by definition sensitive to the 
total mass, and is therefore one of the most promising 
methods to measure the mass profile. Some analyses have 
combined strong- and weak-lensin g or velocity dispersion 
constraints for individual clusters [2 7l [H, [2§, I30II31L [32j 
to derive the profile, and in some cases concentrations 
have been found above the predictions from simulations. 
However, strong lensing is affected by the mass distri- 
bution in the very inner parts of the cluster, and both 
cllipticity of matter and stars from the central galaxy 
have a significant effect on the strong lensing signatures, 
so that these analyses are not necess arily measuring the 
primordial dark matter halo profile [331 ]. Furthermore, 
both strong lensing and X-ray analyses are susceptible 
to selection bias effects when changing the DM concen- 
tration at fixed mass (26j |. The problems in interpreting 
the observations on small scales, where baryons play an 
important role, suggest that we should focus on larger 
scales if we want to compare observations to theoretical 
predictions from N-body simulations. Weak lensing is ar- 
guably the most promising tool that can be used to mea- 
sure the profile out to scales of several h~ 1 Mpc. It has 
the advantage that outside the central region, the dark 
matter distribution is likely to be unaffected by baryons, 
so a comparison against N-body simulations should be 
more reliable. 

Many previous weak lensing analyses have focused on 
individual clusters (for example, [H], HH). Measuring the 
matter distribution of individual clusters has its advan- 
tages, since it allows a comparison with the light and gas 
distributions on an individual basis, and so can constrain 
models that relate the two, such as MOND versus CDM 
[36|. However, since lensing measures a projected surface 
density with a window that extends hundreds of mega- 
parsecs away, other mass perturbation along the line of 
sight will also produce a lensing signal and thus act as a 
source of noise when extracting the cluster density pro- 
file. Some of these structures may be correlated with the 
cluster itself, for example those that are falling into the 
cluster along the filaments connected to the cluster [37l ]. 
and can be defined as part of the cluster profile, while 
other structures may be completely unrelated mass con- 
centrations tens or hundreds of megaparsecs away [38j . 

The measurement of the dark matter profile can there- 
fore be quite noisy for individual clusters. Stacking the 
signal from many clusters can ameliorate this problem, 
since only the mass density correlated with the cluster 
will produce a signal. This way, the measurement deter- 
mines the true average cluster mass profile in the same 
way as defined in simulations. Such a statistical approach 
is thus advantageous if one is to compare the observations 
to theoretical predictions, which also average over a large 
number of halos in simulations. In fact, simulations show 
a significant scatter in the shapes of individual halo pro- 
files [2j, |39( , so stacking many halos will reduce the fluc- 



tuations due to noise caused by uncorrelated structures 
along the line of sight, due to shape measurement noise, 
and due to the shape variations of individual halos. A 
final advantage of stacking is that it allows for the lens- 
ing measurement of lower-mass halos, where individual 
detection is impossible due to their lower shears relative 
to clusters. Individual high signal-to-noise cluster ob- 
servations and those based on stacked analysis of many 
clusters are thus complementary to each other at the high 
mass end, with the stacked analysis drastically increasing 
the available baseline in mass. 

The statistical approach based on stacked clusters has 
been applied to a small number of clusters before (40l.l4l|. 
The Sloan Digital Sky Survey (SDSS) provides an ideal 
dataset for such analysis: it covers a significant fraction 
of the sky containing ~ 10 4 clusters up to z ~ 0.3, pro- 
viding a large volume of the universe over which the clus- 
ters can be observed. The SDSS spectroscopy and multi- 
color imaging enables precise redshift determination, so 
we can determine the profile as a function of true trans- 
verse separation rather than angle. In a previous analy- 
sis, we have used a sample of 43 335 groups and clusters 
as traced by luminous red galaxies (LRGs) [42J to derive 
average mass density profiles of groups with masses be- 
tween 3 x 10 13 /i _1 M Q and 1.3 x 1O 14 /i _1 M . We have 
also performed a halo model analysis of the lensing signal 
of isolated ~ elliptical galaxies in [43[ to ensure that 
they are all in the field, and will augment that sample 
with isolated ~ L* spiral galaxies. In this paper we ex- 
tend these previous analyses to the new sample of 13 823 
MaxBCG clusters presented in [Hj], which extends the 
mass range to ~ 6 x 10 14 H^Mq. We also compare our 
results against an independent analysis of these clusters 
in [45(1 , though both analyses have included objects other 
than the public maxBCG catalog selected in different 
ways. We then combine the elliptical, spiral, LRG and 
MaxBCG analyses to obtain information about the halo 
density profile over a wide range of masses, and compare 
them to theoretical predictions from N-body simulations. 

We begin in section [H] by presenting the theory be- 
hind our measurement, the data, and the analysis method 
used. The results are presented in section IIII1 including 
a discussion of how to compare them with other obser- 
vations and with theory. Our conclusions derived from 
this analysis and comparison with theory are given in 
section [TV] 



II. DATA AND ANALYSIS 
A. Theory 

We follow the same methodology as in [42], so we re- 
fer the reader to that paper for a more detailed descrip- 
tion of the analysis process. In brief, cluster-galaxy and 
galaxy-galaxy weak lensing provide a simple way to probe 
the connection between clusters (or galaxies) and matter 
via their cross-correlation functions £, c i,m(r) (or £ ff , m (r)), 
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which can be related to the projected surface density 



E(R)=p 



1 + £c/,- 



x 



dx (1) 



where R is the transverse separation and x the radial 
direction over which we are projecting. We are ignor- 
ing the effects from the radial window, which is hun- 
dreds of megaparsecs broad and not relevant at cluster 
scales. The surface density is then related to the observ- 
able quantity for lensing, the differential surface density, 



AE(i?) = 7t (fl)£ c - £(< R) - E(JJ), 



(2) 



where the second relation is true only in the weak lens- 
ing limit, for a matter distribution that is axisymmetric 
along the line of sight (which is naturally achieved by our 
procedure of stacking thousands of clusters to determine 
their average lensing signal). This observable quantity 
can be expressed as the product of two factors, a tangen- 
tial shear j t and a geometric factor 



£ c 



D, 
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where Dl and Ds are angular diameter distances to the 
lens and source, and D^s is the angular diameter dis- 
tance between the lens and source. Unless otherwise 
noted, all computations assume a flat ACDM universe 
with TO = 0.27 and J7a = 0.73. Distances quoted for 
transverse lens-source separation are comoving (rather 
than physical) /i _1 kpc, where Ho = 100 h kms" 1 Mpc -1 . 
Likewise, the differential surface density AE is computed 
in comoving coordinates, and the factor of (1 + zl)~ 2 
arises due to our use of comoving coordinates. 

For this paper, we are primarily interested in the 
contribution to the cluster-mass or galaxy-mass cross- 
correlation from the halo profile itself, rather than from 
neighboring halos (halo-halo term), and hence 
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The halo-halo term for clusters in host halos can be 
modeled simply using the cluster- or galaxy-dark mat- 
ter cross-power spectrum as in, e.g., [46[, and is only 
important for R >~ 2 /i _1 Mpc. Nevertheless, we com- 
pute this component and include it in the model as a 
fixed term which we obtain by computing first the mass 
of the clusters, deriving the corresponding halo bias us- 
ing the bias-mass relation EH and using the linear 
power spectrum multiplied with the halo bias to obtain 
the halo-halo term. This procedure has been shown to 
work well in comparison to simulations [46[ , so we use the 
bias-mass relation from that paper but with a$ = 0.75 
and Q m = 0.27, computed at the mean redshift of each 
sample. 

We can model the one-halo term for each sample as a 
sum of the stellar component, only important on scales 



below 
El: 



lOO/i 1 kpc, and an NFW dark matter profile 



p(r) 



(r/r a ) (1 -I- r/r s 
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It is convenient to reparametrize it by two parameters, 
concentration c = r V i r /r s and virial mass M. The virial 
radius r V i r and p s can be related to M via consistency 
relations. The first is that the virial radius is that within 
which the average density is equal to 200 times the mean 
density: 
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The second relation, used to determine p s from M and c, 
is simply that the volume integral of the density profile to 
the virial radius must equal the virial mass. The NFW 
concentration c is a weakly decreasing function of halo 
mass, with a typical dependence as 
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with (3 ~ 0.1 @, H^, |49j], making this profile a one- 
parameter family of profiles. The normalization depends 
on the nonlinear mass, which is the mass within spheres 
in which the rms fluctuation in the linear regime is 1.68; 
for the typical range of cosmological models, one expects 
c w 5 - 6 at Mo = 1O u /i _1 M0. 

As an example of how the lensing signal varies with 
concentration, in Fig. Q] we show the predicted NFW 
lensing signal as a function of transverse separation for 
M 20 o6 = 10 12 (galaxy scale) and 10 u /i -1 Mg (cluster 
scale) halos, for several plausible concentration values. 
The vertical line shows the minimum scale used for our 
fits, and it is clear that the lensing signal above those 
scales can differentiate between different concentration 
values. This discriminating power stems in part from the 
fact that the lensing signal reflects the differential surface 
density, which draws information from smaller scales to 
larger scales. 

We fit the data to the model assuming a spherical NFW 
profile. In [46| . it was shown that spherical NFW pro- 
files do an excellent job at describing the stacked lensing 
signal from simulations for a variety of masses, with the 
masses and concentrations of the best-fit profiles related 
to the real masses and concentrations in the simulations 
in a particular way (to be discussed further below). The 
fit x 2 values were very good when the errorbars used on 
the simulated lensing signal were one-tenth of our cur- 
rent errorbars, so henceforth we consider only spherical 
NFW fits rather than trying to account for the averaging 
of triaxial halos. 

Although we could include the effects of the central 
galaxy to make the model accurate on scales below 20% 
of the virial radius, there are reasons why we should ex- 
clude this information from the fits if we want to compare 
against N-body simulations. First, the baryonic effect 
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FIG. 1: The NFW profile theoretical lensing signal AS for 
two mass scales. In the upper plot, with M = 10 12 Mq, 
the concentration is varied from 7 to 10 to 13. In the lower 
panel, with M = lO 14 ft _1 Af0, it is varied from 3 to 5 to 7. 
The curves on each panel are normalized to the same virial 
mass M2oot.. Vertical lines show the minimum scale for our 
fits. 



we assume may not be completely accurate because of 
uncertainties in the stellar mass to light ratio, and in 
the dark matter response to the presence of baryons and 
stars forming out of them. The latter is often modeled 
as adiabatic contraction [Hoi . IHlj . but this prescription 
may inaccurately describe the actual effect, depending 
on the formation history of the galaxy, group, or cluster 
in question (52| . The second, possibly more important 
effect is that for cluster samples such as MaxBCG, the 
cluster center cannot always be reliably determined us- 
ing the optical information. This uncertainty has two 
causes: first, that a brightest cluster galaxy (BCG) may 
not be located at the deepest part of the cluster potential, 
but may be offset from it due to (for example) pertur- 
bations from infalling satellites; and second (and more 
importantly), that the maxBCG algorithm may choose 
the wrong BCG. Studies comparing the BCG position to 
the center defined by either X-ray intensity or by average 
satellite velocity have found that the typical displacement 
is about 2-3% of the virial radius when the BCG is prop- 
erly identified @, [H [H[ . The last of these studies 
finds that for about 10% of BCGs, the displacement ex- 
tends to above 10% of the virial radius. It is interesting 
to note that selecting blue core BCGs significantly re- 
duces the displacement to below 1% of the virial radius 
55]. However, while the multi-band SDSS aperture pho- 
tometry could be used to select blue-core BCGs, we do 
not undertake this approach for this analysis because of 
decreased S/N. 



Another study that includes red galaxy photometric 
errors (i.e., both effects rather than just the first) finds 
that the median displacement is 10% of the virial ra- 
dius [5(| . If the assumed center is displaced from the 
true center, then this can have a significant effect on the 
density distribution in the inner parts, mimicking a halo 
profile with a lower concentration [57], [H| , while in the 
outer parts the effect is significantly smaller. To avoid 
this problem, we use a minimum scale for the fits of 500 
/i _1 kpc for the maxBCG clusters (typically half the virial 
radius) , and present evidence that this approach is robust 
to the effects of centroiding problems, at the expense of 
decreased statistical power. 

We note that this problem has also been addressed in 
a statistical manner, using mock catalogs to determine 
both the fraction of clusters affected by centroiding prob- 
lems, and also the distribution of projected offsets from 
the true cluster center for those that are centroided wrong 
PH - Using this correction, one can in principle correct 
for the effect and fit to significantly smaller scales; how- 
ever the correction is quite dependent on the content of 
the mock catalogs, and any deviation from reality may 
invalidate it. For example, the mock catalogs suggest 
that misidentification of the center is more of a prob- 
lem at lower halo masses, while our visual inspection of 
clusters suggests the opposite, in the sense that the most 
massive clusters (with richness A^oo > 80) have signif- 
icant bimodality and substructure, as expected from a 
hierarchical cosmological model where the most massive 
clusters formed the latest. As a result, we exclude these 
most massive clusters from the analysis entirely. For the 
remaining sample, we choose the different approach of 
avoiding the inner parts of the cluster (where the effect 
is strongest), thereby decreasing the systematic error at 
the expense of an increase in statistical error. 



B. Data 

The data used here are obtained from the SDSS 59], 
an ongoing survey to image roughly 7r steradians of the 
sky, and follow up approximatel y o ne million of the de- 
tected objects spectroscopically [60t l6ll . |62| . The imag- 
ing is carried out by drift-scanning the sky in photomet- 
ric conditions [H, 0], i n fi ye bands (uariz) [6p . |66| us- 
ing a specially-designed wide- field camera [63]. These 
imaging data are used to create the source catalog that 
we use in this paper. In addition, objects are targeted 
for spectroscopy using these data [68| and are observed 
with a 640-fiber spectrograph on the same telescope [6!| . 
All of the data are processed by completely automated 
pipelines that detect and measure photometric proper- 
ties of objects, and astrometrically calibrate the data 
0, IzOi. Th e SDSS has had seven major data re- 
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In the subsections that follow, we describe the lens and 
source samples. 
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1. MaxBCG cluster lenses 

Our highest mass lens sample consists of 13 823 
MaxBCG clusters @ H3|, which are identified by the 
concentration of galaxies in color-position space, using 
the well known red galaxy color-redshift relation [8l| . 
The sample is based on 7500 square degree of imaging 
data in SDSS. There is a tight mass- richness relation 
that has been established using dynamical information 
across a broad range of halo mass [82| • This is also con- 
firmed with the lensing analysis, as shown below and in 

EEIlSIll- 

The redshift range of the MaxBCG catalog is 0.1 < z < 
0.3, and the upper cutoff ensures that the lenses still have 
a sufficient number of sources behind them. Within these 
redshift limits, the sample is approximately volume- 
limited with a number density of 3 x 10~ 5 (/i/Mpc) 3 , ex- 
cept for a tendency towards higher number density at the 
lower end of this redshift range [84{ . The main tracers of 
halo mass provided by the MaxBCG team are a rescaled 
richness A200 (number of red galaxies above 0.4L*), to- 
tal luminosity L200 (the luminosity in those red galaxies, 
including the BCG) and BCG luminosity Lbcg- 

In this paper we use richness -/V200 as a primary tracer 
of halo mass. However, we observe (based on 84]) that 
the lower mass end of the public catalog overlaps with the 
LRG samples used for this work, so we eliminate those 
clusters with the lowest richness, A^oo = 10 and 11 (1/3 
of the public catalog), and split the remainder into six 
narrow A^oo bins as shown in Table U (except for those 
with A^oo > 80, as described in Sec. Ill A[) . The widest of 
these spans a factor of 1.4 in A^ocb so if we assume the 
mass is proportional to this observable, then the mass 
bin is a factor of 1.4 wide without any mass-observable 
scatter. In reality such scatter exists, but even a factor of 
two scatter in the mass at fixed A^oo would give a mass 
distribution less than an order of magnitude wide. As we 
have shown in [461 ] , when fitting for a single halo mass on 
a stacked sample that is less than an order of magnitude 
wide, the best-fit mass is a good proxy (within ~ 10%) 
for the true mean mass of the sample. 



2. Spectroscopic LRGs 

For the next sample with lower average halo mass, we 
use the spectroscopic Luminous Red Galaxy (LRG) sam- 
ple [6(3], including area beyond Data Release 4 (DR4). 
The total area coverage for this spectroscopic sample 
is 5154 square degrees, as available in the NYU Value 
Added Galaxy Catalog (VAGC, [IH) at the time of the 
original publication of our lensing work with this sample 

For the LRG sample, we define luminosities using r- 
band model magnitudes, extinction-corrected using red- 
dening maps from [86| . We apply a k+e-correction (com- 
bined /c-correction and correction for evolution of the 
spectrum) to redshift zero using stellar population syn- 



thesis code from [87j ; this magnitude is denoted Afo.o,. 
to distinguish it from the magnitudes used for the Main 
sample, which are corrected to z = 0.1. In our origi- 
nal work, the LRGs were split into two luminosity bins, 
Mo.o r > —22.3 and Mo.o r < —22.3. To reduce the over- 
lap between the LRGs and the cluster samples to < 5%, 
we use the fainter bin in its entirety but only use the 
—22.6 < Ala.a r < —22.3 subset of the brighter sample. 

These LRG samples were derived from the full LRG 
sample after eliminating 15% of the sample using a 
cylindrical density estimator, designed to avoid satellite 
galaxies due to their extra lensing signal from the host 
halo. Specifically, the LRGs were each required to be 
the only or the brightest LRG in a cylinder with radius 
R = 2/i~ 1 Mpc and line-of-sight separation 5v — ±1200 
km s . This cut is conservative, in the sense that for 
typical groups and low-mass clusters, we may have ex- 
cluded some host galaxies; however, the host sample pu- 
rity is sufficiently important for this analysis that we tend 
towards the conservative side. More information about 
the LRG samples is available in Table [TJ 



3. L, lenses 

Finally, we include lower luminosity samples from [43| 
that have been shown with a robust environment esti- 
mator to consist of field galaxies. For those samples, we 
used galactic extinction-corrected r-band Petrosian mag- 
nitudes, fc-corrected to z = 0.1 using KCORRECT v4_1_4 
[H], denoted Mo.i r . 

The samples used here correspond to L3, L4, and 
L5faint isolated ellipticals from [43|], where the isolated 
ellipticals are the half of the elliptical sample at those lu- 
minosities determined to be in the field using a cylindrical 
density estimator. We also include a sample not shown 
there but from the same data, of L3 and L4 isolated spi- 
rals, or 85% of the spirals in those luminosity bins. For 
reference, L* is within L4, so our samples range from 
slightly below to slightly above L*. In this context, the 
elliptical and spiral samples are chosen using the SDSS 
frac_deV parameter, which determines whether the light 
profile is closer to a de Vaucouleurs or exponential pro- 
file. More details of these five lens samples are shown in 
Table D 

We do not use the fainter samples (LI and L2) from 
previous work because the detection significance of the 
weak lensing signal is low, so they cannot constrain the 
halo concentration. We also avoid the brighter ellipti- 
cal samples because they overlap with the spectroscopic 
LRGs, and the brighter spiral samples because they are 
nearly empty. 

4- Lensing sources 

The source sample used is the same as that originally 
described in [8!|. This source sample includes over 30 
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TABLE I: Summary of all lens samples used in this paper. Note that observable definitions (such as luminosity) are not 
necessarily the same for each parent sample; see the relevant subsection for details. Additional sample cuts are indicated by 
the parent samples, which typically are defined using flux and color cuts. 
Parent sample Isolation method Observable cut Redshift cut (z) N Name 



MaxBCG clusters MaxBCG method 



12 < A200 < 13 
14 < iV 2 oo < 19 
20 < A200 < 28 
29 < A200 < 39 
40 < A200 < 54 
55 < A200 < 79 



0.1 < z < 0.3 0.22 



2531 
3372 
1618 
614 
248 
109 



Spectroscopic LRGs Cylindrical 



Mo.o r > -22.3 
-22.3 > Mo.o„ > -22.' 



0.15 < z < 0.35 0.24 27700 
10536 



MAIN spectroscopic Cylindrical 



-20 < Mo.i r < -19, frac_deV> 0.5 
-21 < Mo.i r < -20, frac_deV> 0.5 
-21.5 < Mo.i r < -21, frac_deV> 0.5 
-20 < Mo.i r < -19, frac_deV< 0.5 
-21 < Mo.i r < -20, frac_deV< 0.5 



2 > 0.02 0.07 20150 L3 elliptical 
0.10 46130 L4 elliptical 
0.13 23485 L5faint elliptical 
0.07 38640 L3 spiral 
0.10 42235 L4 spiral 



million galaxies from the SDSS imaging data with r-band 
model magnitude brighter than 21.8, with shape mea- 
surements obtained using the REGLENS pipeline, in- 
cluding PSF correction done via re-Gaussianization (90j 
and with cuts designed to avoid various shear calibra- 
tion biases. In addition, there are also uncertainties due 
to photometric redshifts and/or redshift distributions of 
background galaxies, which were originally calibrated us- 
ing DEEP2 Groth strip data, as well as due to other is- 
sues affecting the calibration of the lensing signal, such 
as the sky subtraction uncertainties, intrinsic alignments, 
magnification bias, star-galaxy separation, and seeing- 
dependent systematics. The overall la calibration uncer- 
tainty was estimated to be eight per cent [8fj| , though the 
redshift calibration component of this systematic error 
budget has recently been significantly decreased due to 
the availability of more spectroscopic data [9l[ . The cal- 
ibration mainly affects the mass estimation rather than 
the derived density profiles, so it is not of significant con- 
cern for this paper due to the weak dependence of con- 
centration on mass. 

An additional concern is the relative calibration of 
lensing measurements from different lens samples that 
were published in different papers. If there are calibra- 
tion differences between these measurements, then the 
power-law scaling of the c{M) relation might be mises- 
timated. However, it seems unlikely that there can be 
significant calibration differences between the different 
measurements, for two reasons. First, the same version 
of the source catalog was used for each one. This suggests 
that any shear calibration or star/galaxy separation is- 
sues are the same for each measurement. There may be 
very slight variation due to the different mean lens red- 
shifts, which changes the effective mean redshift of the 
sources, but our previous tests of the source catalog for 
relative shear calibration as a function of apparent mag- 
nitude and size rule out changes in the calibration of 
the shear that are significant relative to the la statistical 
error on these measurements. Second, we have rigorously 



tested the calibration of the source redshift distributions 
using the zCOSMOS and DEEP2 spectroscopic samples, 
and found that the calibrations for the lens redshift dis- 
tributions used here are the same within several percent, 
which is again smaller than the la measurement error on 
the lensing signals used here In short, while there 

is some small uncertainty (discussed above) in the ab- 
solute lensing signal calibration, we have little reason to 
believe there is any significant discrepancy between the 
calibrations for the different lens subsamples. 



C. Analysis 

Calculation of the signal is described in detail in [42| . 
Briefly, we compute the weights based on noise and 
redshift information for each lens-source pair, summing 
them using a minimal variance estimator. We compute 
the signal around random points and subtract it from 
the signal around real lenses to eliminate contributions 
from systematic shear. The signal must be boosted, i.e. 
multiplied by B(R) = n(R)/n ran d(R), the ratio of the 
number density of sources relative to the number around 
random points, in order to account for dilution by sources 
that are physically associated with lenses, and therefore 
not lensed. The former correction is only important on 
scales above those used in this paper (> 5 /i _1 Mpc) and 
the latter on scales below 20% of the virial radius, which 
we do not use for the fits in this paper. To determine er- 
rors on the lensing signal, we divide the survey area into 
200 bootstrap subregions, and generate 2500 bootstrap- 
resampled datasets. We note that the effects of non-weak 
shear, magnification bias, sky subtraction and intrinsic 
alignments, discussed in more detail in [42l. [89j. are neg- 
ligible on the scales used in this paper. 

The errors determined from the bootstrap are used in 
plots of the signal with errors; however, the bootstrap co- 
variance matrices can be quite noisy and therefore inap- 
propriate to use for weighting in the x 2 minimization fit- 
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ting. To avoid this problem, we determine analytic, diag- 
onal covariance matrices (including shape noise), which 
we have shown in [8§| to be a less noisy version of the 
bootstrap covariance matrices with agreement in size of 
the errors at the ~ 10% level. These covariance matri- 
ces, which are far less noisy, are used to perform the fits 
on each bootstrap-resampled dataset. The distributions 
of output parameters from all the bootstrap-resampled 
datasets are used to determine errors on the fit parame- 
ters. 

The fits are for eleven parameters, using the lensing 
signal from the thirteen lens samples in Table HI 

• The normalization and slope of the concentration- 
mass relation, Eq. [7] (2 parameters). 

• The normalization and slope of the relation be- 
tween mass M2006 and MaxBCG richness A^oo (2 
parameters): 

M 2006 = M (^y (8) 

• The masses of the two LRG and five lower luminos- 
ity samples (7 parameters). 

As shown in Eq. the concentration is expected to 
scale with redshift. Since the samples are at different 
mean redshifts, we use the expected redshift scaling to 
fit for a normalization cq(z = 0.22) (the mean, lensing- 
weighted redshift of the maxBCG sample) . 

III. RESULTS: LENSING SIGNAL AND FITS 
TO NFW PROFILE 

We perform the analysis described in section III CI on 
the lensing signal AS for all the lens samples. Fig. [2] 
shows the lensing signal AS with the best-fit model in the 
joint fits for c 2 oob(M 2 oob, z = 0.22), for Af 2 oofc(A r 2 oo) for 
the maxBCG clusters, and best-fit masses for the spec- 
troscopic LRGs and the lower luminosity galaxy samples. 
Fig. [3] shows the signal and best-fit model for the lower 
luminosity samples, and Fig. [4] shows the same for the 
two LRG samples. We see that for the maxBCG sample 
there is a strong lensing signal over the entire range of 
richness, and the lensing signal is increasing with richness 
as expected. The best-fit parameters for concentration 
and mass (to be described below) clearly provide a good 
fit to the data in all lens samples, as will also be evident 
in the % 2 values for the fits, discussed below. 

A. Description of fits 

Table [TT| shows the best-fit parameters for the power- 
law relations c(M) and maxBCG M(N2oo), using sev- 
eral approaches to the fits as described below. The ta- 
ble includes information about the fit minimum radii and 



whether offsets of BCGs were accounted for in the fitting. 
It also shows the \ 2 P er degree of freedom, indicating 
that all fits shown provide a reasonable fit to the data. 
As will be discussed in detail later, Fit 2 in that table is 
used for any plots that include a best-fit model, because 
the fit procedure represents the best possible trade-off 
between systematic error due to cluster centroiding er- 
rors and statistical error. Table IIIII shows the best-fit 
masses for the two LRG lens samples and the five lower 
luminosity samples for fit 2 in Table HH 

The tables include the results for several types of fits 
with varying minimum radii for the maxBCG sample; 
these results were shown to test several possible issues. 
In all cases, the maximum fit radius for the maxBCG 
clusters was 3.0 /i _1 Mpc; the fits for the LRGs used 0.1- 
2 /i _1 Mpc; and the fits for the lower luminosity samples 
used 0.04-0.5 /i -1 Mpc. The maximum fit radius was cho- 
sen so that the halo-halo term remains small compared to 
the one halo term. When computing the halo-halo term 
we assume the cosmological model and redshift of these 
samples described in Section III Al For example, for the 
lower mass LRG sample, a mass of 3 x 10 13 /i _1 M© gives 
a bias of 1.8. We compare this number against the results 
in [92j, which gives a clustering amplitude for this sam- 
ple relative to L* of 1.85. The highest richness maxBCG 
sample, approaching M ~ 6 x 1O 14 /i -1 M0, has a bias of 
5.5 in this model. 

We justify our neglect of the stellar term by consider- 
ing that a point mass at the origin gives a lensing sig- 
nal AS = M point / (irR 2 ). For the L* samples, typical 
stellar masses are a few xl0 lo /i -1 M© and the minimum 
fit radius is 0.04 ft, _1 Mpc, which implies that the signal 
from the stellar component at that radius is subdominant 
compared to the observed signal of tens of HMq /pc 2 . For 
the LRGs, the stellar mass is higher by factors of a few, 
but the minimum fit radius is 2.5 times higher, which 
again makes the stellar component subdominant on the 
scales used for the fit. Finally, for the maxBCG sample, 
given point masses of order 10 12 /i _1 M Q and at least 0.2 
/i _1 Mpc for the minimum fit radius, it is difficult to ar- 
range for the stellar term to be more than 5 per cent of 
the predicted signal from the NFW dark matter halo. 



1. Basic fits 

We begin by discussing the basic fits 1-3, which use 
minimum fit radii of 0.2, 0.5, and 1 ft, Mpc, with- 
out explicitly accounting for failure to properly centroid 
BCGs. We first look for systematic effects due to BCG 
offsets by comparing the fit results with different min- 
imum radii. The signature of BCG offsets would be 
a lower concentration for a lower minimum fit radius; 
we see no definitive sign of any such problem on the 
scales used here. For the case with minimum fit ra- 
dius of 0.5 /i _1 Mpc, which we adopt as a compromise 
between minimizing the systematic error error due to 
this effect while maximizing the statistical contraining 
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FIG. 2: The lensing signal AE and model prediction for the best-fit model parameters (fit 2 in Table Hljl . for the maxBCG 
sample split into 6 richness bins. 



power, the concentration at Af2oo6 = 1O 14 /i -1 M0 is 
co = 4.6 ± 0.7, with a power-law scaling with mass of 
—j3 — — 0.13±0.07. The mass normalization at A^oo = 20 
is M 20 o6 = (1.56 ±0.12) x lO^^Mg, with a scaling 
with richness of 1.15 ± 0.14. The \ 2 of 266 for 282 de- 
grees of freedom indicates that the fit is acceptable, as 
shown in Figs.(2HH Our minimum maxBCG richness bin, 
at (-/V200) = 12.4, therefore has mean mass and concen- 
tration in this model of M 2 oob = 0.9 x l0 14 h~ 1 M Q and 
C200h = 4.7. 

The best-fit masses for the LRG and lower luminos- 
ity samples are shown in Table IIII1 and as we antici- 
pated, the LRGs are in a lower mass range that does not 
overlap that of the maxBCG sample. The spectroscopic 
LRGs with Mo.o r > —22.6 therefore predominantly trace 



group-scale halos, below the cluster-scale masses of the 
maxBCG sample with -/V200 > 12. We also confirm our 
previous results [431 ] that, at fixed r band luminosity, iso- 
lated L* spirals have a lower mass than isolated L* el- 
lipticals. The signal detections are higher significance 
than the gm/M values would indicate for L3: the error 
distributions determined using the bootstrap-resampled 
datasets are non-Gaussian and well-separated from zero. 



2. Fits with offsets 

We next consider the other sets of fits in Table HH Fits 
4-6 differ from 1-3 only in the inclusion of a prescrip- 
tion given in [45j to correct for the centroiding problem 
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FIG. 3: The lensing signal AE and model prediction for the best-fit model parameters (fit 2 in Table ITT)) . for the field low 
luminosity sample split into 5 morphology and luminosity bins (left column: ellipticals; right column: spirals). 



in the maxBCG catalog (with no change in the way the 
lower mass samples were handled) . This prescription, de- 
rived from mock catalogs, is described in detail in [45j . 
In brief, it has a richness-dependent fraction of misiden- 
tified BCGs (from 30% at low richness to 20% at high 
richness), and those that are misidentified have a Gaus- 
sian distribution of projected separation from the true 
centroid, with a scale radius of 0.42 h~ 1 Mpc. 

The key point to consider is that if the model used 
to account for centroiding problems is correct, then the 
best-fit concentrations should be independent of the min- 
imum scale used. We see that this is not the case: the 
concentrations are significantly elevated compared to pre- 
vious results for minimum fit radii of 0.2 and 0.5 ft, Mpc, 
and differ from the results with a minimum fit radius of 



1 ft, _1 Mpc, for which (as one expects) the offsetting pre- 
scription does not significantly change the results given 
that the minimum fit radius is well outside the offset 
scale. Thus, we suggest that this particular offset pro- 
cedure, which is determined using mock catalogs in [45j, 
may in fact overcompensate for the true level of BCG 
offsets. This would not be too surprising, given that it 
depends sensitively on the contents of the mock catalogs, 
and is only applicable if they are a very true representa- 
tion of the real world. At the high-mass end, the com- 
parison of maxBCG versus X-ray centroids in [44j can be 
used to evaluate the offset procedure derived from mock 
catalogs. For the derived 20 per cent failure rate, with a 
Gaussian scale length of 0.42 /i _1 Mpc for the projected 
offset, we would expect 6.9 and 8.7 per cent of the clus- 
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FIG. 4: The lensing signal AE and model prediction for the 
best-fit model parameters, for the host LRG sample in two 
luminosity bins as labeled on the plot. 



ters to have centroids with projected offsets of 0.25-0.5 
and 0.5-1 /i _1 Mpc, respectively. Out of the 87 clusters 
with X-ray matches within 1 h~ 1 Mpc, figure 14 in [44| 
shows 11 (12.6 per cent) and 4 (4.6 per cent) in those 
two ranges of projected offsets. While these numbers are 
not formally inconsistent with the expectations from the 
mocks (e.g., for the 0.5-1 /i _1 Mpc offset range, the one- 
tailed P(n < 4) = 0.12), had they been used to derive an 
empirical model for the offset, it would entail a smaller 
correction than that derived from the mocks at the scales 
which are most relevant for our adopted fit 2. Since there 
may be uncertainty in the X-ray centroids for disturbed 
or merging clusters, even the lower estimate in [HJ may 
be an overestimate. 

The fact that [56| find similar levels of offsets for pho- 
tometric LRGs from the centroids of cluster X-ray dis- 
tributions does not contradict this conclusion: an exam- 
ination of those results suggests that they only apply for 
masses above ~ 5 x 10 14 /i -1 M Q (after converting mass 
definitions) , which essentially corresponds to the top two 
maxBCG bins considered here. Furthermore, the mean 
redshift for that study is higher than here, which implies 
different levels of photometric redshift errors and there- 
fore different levels of BCG misidentification. Thus, it 
may not be applicable in detail here even for the top two 
bins. Finally, the use of higher mass clusters at higher 
redshifts should give a higher fraction of clusters that 
experienced recent mergers, so again we expect an over- 
estimate of the effect for the full maxBCG sample. The 
true level of offsets in the four lower bins, which dom- 
inate the c(M) fits, is poorly constrained from the real 
data. 



3. Fits with fixed power-law scaling 

Next, to investigate the correlation between cq and 
(3, we consider fits 7-12, which are the same as 1-6 ex- 
cept with a fixed f3 = 0.1 (the theoretical value). While 
Table ILT1 gives the formal correlation coefficient between 
these parameters from the fit covariance matrix, it is in- 
structive to explicitly fix f3 = 0.1 to see the effect on cq. 
As shown, there is some degeneracy between cq and (3, 
but a comparison of e.g. fits 2 and 8 suggests that this 
degeneracy is not very strong for the fit we have selected 
as our main result. Note that for fits 10-12, fixing (3 has 
somewhat ameliorated the discrepancy between the fits 
using the offsetting procedure with different minimum 
scales. Nonetheless, the trend towards increasing c with 
decreasing minimum scale is suggestive of possible over- 
compensation for the true level of the problem in reality. 

As stated previously, we choose fit 2 as our most ro- 
bust result, but the difference between this and the other 
fits suggests a systematic error that is comparable in size 
to the statistical error. Without better knowledge of the 
true level of offsets of the maxBCG from the cluster cen- 
troids (in the statistical sense) , it is impossible to reduce 
this systematic error. However, for masses below a few 
xl0 14 /i _1 M Q , where X-ray cluster data are difficult to 
obtain with sufficient resolution at these redshifts, there 
is no clear, simple way to observationally constrain these 
offsets at this time. 

4- Systematic tests of fitting procedure 

Next, we discuss a few fits not included in the table 
that were designed to test for systematic errors. [§3[ 
find, using the Millenium simulation, that DM halos can 
be more properly described using the Einasto profile than 
the NFW profile, 

/0 (r) = Ps e(- 2 /«)[(^)"-4, ( 9 ) 

where a has a weak mass dependence with a value around 
0.15. Thus, fits to NFW profiles in simulations can lead 
to different concentrations depending on the scale used 
for the fits. While the differences between Einasto and 
NFW profiles are most significant well within the scale 
radius, where we do not probe using weak lensing, we 
nonetheless test that our results are insensitive to the 
choice of Einasto versus NFW profiles. When doing a 
fit that is comparable to fit 2 but with NFW profiles 
replaced by Einasto profiles, we find the best-fit masses 
to be preserved, and the best-fit c = 4.5 ± 0.7 and (3 
0.12 ± 0.7. The changes from fit 2 are well within the 
la error, so we conclude that possible errors in best-fit 
parameters due to the use of NFW rather than Einasto 
profiles are insignificant. 

We also perform fit 2 using three broader bins in ./V200 
instead of the six narrow bins used for the main re- 
sults. We find that the recovered c(M) relation is vir- 
tually unchanged; the main parameter that varies is 7 in 
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the M(7V2oo) relation, which becomes shallower by 1.5er. 
This result is as expected from [46J: the broad bins are 
most problematic at the highest mass end, where they 
reduce the best-fit mass. Since we are not fitting the 
bins for masses individually, but rather are fitting for a 
power-law relation, the exponent of this relation is con- 
sequently reduced. Given the size of this shift, and the 
fact that our narrower bins used for the main analysis 
should contain mass distributions less than an order of 
magnitude wide, we do not ascribe significant systematic 
error in the concentration-mass relation to the default 
bin size. 

Finally, because of potential centroiding systematics in 
the maxBCG lensing sample that should not be present 
for the lower luminosity or LRG samples, we performed 
the fits without the maxBCG samples entirely. In that 
case, we find that is quite poorly constrained, so we fix 
it to 0.1 and compare against fit 8 in Table HH which also 
has (3 = 1 and only differs in that it includes the maxBCG 
sample. In this case, we find the best-fit Co = 4.7 ± 0.7, 
entirely consistent with the results in the table. This 
result suggests that our choice of minimum fit radius has 
minimized systematic error due to maxBCG centroiding 
errors to be well within the statistical error. 



B. Concentration-mass relation 

In Fig. [5l we show the best-fit c(M) relation from 
fit 2, with a lcr error region defined by the fits to 
fifty bootstrap-resampled datasets. As shown, the 
concentration-mass relation is best constrained from 
lO 13 -lO 14 /i _1 Af , due to the interplay between higher 
mass increasing the lensing signal versus higher mass 
meaning a lower number density (and therefore higher 
measurement error). We emphasize that this is the c(M) 
relation at z ~ 0.22, so in the simplest approximation of 
no mergers, the normalization at z — should be higher 
by about 20 per cent. As shown, the range dominated by 
the ~ L» samples (10 12 /i _1 M Q ) yields a concentration of 
10 ±3, as expected theoretically. At 6 x 10 14 /i _1 A/ Q , the 
top end of the maxBCG sample, the constraint is 4±1. A 
constant concentration-mass relation is just barely per- 
mitted at the 2cr level. The red points on the plot are the 
best-fit concentrations and masses for the individual lens 
samples when we fit for c and M for each one without re- 
quiring a power-law c(M) relation. The consistency with 
this power-law indicates that within the errorbars, the 
c(Af ) power-law is indeed a good fit to the data. 



C. Comparison against previous observations 

We can compare these results to our previous lensing 
results based solely on the LRG sample [42]. In that 
case we found C2oofc = 5.2 ± 0.6 at the pivot mass of 
~5x 1O 13 /i _1 A/ , with weak constraints on the slope of 
the mass concentration relation given the narrow mass 




FIG. 5: The best-fit c(M) relation at z = 0.22 with the lcr 
allowed region indicated. The red points with errorbars show 
the best-fit masses and concentrations for each bin when we 
fit them individually, without requiring a power-law c (M) re- 
lation. The blue dotted lines show the predictions of [3y] for 
our mass definition and redshift, for the WMAPl (higher) and 
WMAP3 (lower) cosmologies. The prediction for the WMAP5 
cosmology falls in between the two and is not shown here. 



range traced by LRG halos. This number is in good 
agreement with our fiducial value C2oob = 4.7 ± 0.7 at 
10 14 /i _1 M©, which gets increased by 10 per cent when 
going to the lower LRG masses. The LRG sample is one 
of the three samples used here and we follow essentially 
the same analysis, so the agreement is to some extent 
expected. 

Next we compare our results against the weak lensing 
determination of c(Af) in [45| . which differs from ours in 
several notable points: (1) we use the maxBCG sample 
to cover the range of masses from 0.8 to 6 x 10 li h~ 1 M@, 
whereas they use a proprietary version of the catalog that 
extends roughly 1.5 orders of magnitude lower in mass; 
(2) we include several additional mass tracers extending 
the halo mass range a factor of ~ 10 lower than in [45] 
with very different selection criteria; (3) we avoid scales 
that are affected significantly by BCG centroiding prob- 
lems, rather than using a correction procedure derived 
from mock catalogs; and (4) the photometric redshifts 
that they use to determine source redshifts and therefore 
normalize the lensin g si gnal suppress the lensing signal 
by ~ 15-20 per cent |9lj. 

In their table 10, they show fit results for power-law 
relations between mass and richness, and concentration 
and richness. We consider their result for Miso&> which 
should differ from our results with A/200& by only a few 
per cent. At iV 2 oo = 20, they find a best-fit mass of 
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TABLE II: Results of fits to lensing signal for several different minimum fit radii for the maxBCG sample split into six richness 
bins, with and without the prescription to account for BCG offsets described in section III Al Fits 1-6 have j3 (the power-law 
slope for the c(M) relation) free; it is fixed to the theoretical value of 0.1 in fits 7-12. 



Fit 


Rmin 


Offsets 


CO 





Corr(c ,/?) 


Mo 


7 y?/v 




/i^Mpc 














1 


0.2 


No 


4.2 ±0.5 


0.16 ±0.07 


-0.37 


1.56 ± 0.12 


1.20 ±0.14 314.8/336 


2 


0.5 


No 


4.6 ±0.7 0.13 ±0.07 


-0.63 


1.56 ± 0.12 


1.15 ± 0.14 266.2/282 


3 


1.0 


No 


4.0 ±0.9 


0.18 ±0.09 


-0.75 


1.53 ± 0.12 


1.23 ± 0.18 215.3/240 


4 


0.2 


Yes 


5.8 ±0.5 


0.03 ±0.07 


-0.44 


1.69 ±0.13 


1.07 ±0.15 317.1/336 


5 


0.5 


Yes 


5.8 ± 0.7 0.03 ±0.07 


-0.67 


1.69 ± 0.12 


1.08 ±0.14 270.4/282 


6 


1.0 


Yes 


4.2 ± 0.9 


0.16 ±0.09 


-0.75 


1.61 ± 0.12 


1.22 ±0.18 214.3/240 


7 


0.2 


No 


4.3 ±0.4 


0.1 




1.53 ±0.11 


1.17 ±0.13 315.9/337 


8 


0.5 


No 


4.8 ±0.6 


0.1 




1.55 ±0.11 


1.14 ±0.13 266.2/283 


9 


1.0 


No 


4.6 ±0.8 


0.1 




1.50 ±0.11 


1.19 ±0.17 215.7/241 


10 


0.2 


Yes 


5.5 ±0.4 


0.1 




1.73 ±0.12 


1.10 ±0.14 318.2/337 


11 


0.5 


Yes 


5.3 ±0.6 


0.1 




1.72 ±0.11 


1.11 ±0.13 271.3/283 


12 


1.0 


Yes 


4.7 ±0.8 


0.1 




1.58 ± 0.11 


1.19 ±0.17 214.5/241 



TABLE III: Best-fit masses for LRGs and lower luminosity 
lens samples from the simultaneous fits to the lensing signal 
corresponding to fit 2 in Table ITT1 



Sample 


Mass M 200 i, (lO iz ft _1 M ) 


LRG (fainter) 


30 ±3 


LRG (brighter) 


56 ±5 


Elliptical L3 


1.0 ±0.4 


Elliptical L4 


1.8 ±0.4 


Elliptical L5faint 


4.5 ±0.8 


Spiral L3 


0.9 ±0.3 


Spiral L4 


0.8 ±0.3 



1.2 x 10 14 h~ 1 Al Q . This number is lower than our result 
in Table ITTI bv 22 per cent; however, there is a straight- 
forward reason for this difference. Using a large spectro- 
scopic training sample, we have recently shown that the 
photometric redshift algorithm used for sources in that 
work leads to a suppression of the lensing signal for these 
lens redshifts of ~ 20% [9l( and hence leads to about 30% 
underestimation of mass. They also find a steeper scaling 
with richness, 1.3 instead of 1.15 as in our work. This 
difference can be explained by the fact that they use a 
much larger range of richnesses for the fit, A^oo > 1 in- 
stead of > 12 as in our work. It is apparent from their 
figure 11 that if one restricts to A^oo > 12, the best-fit 
power-law should be shallower than the result for their 
full richness range, roughly consistent with our result. 

We also compare against their results for concentra- 
tion as a function of richness. They find cisob = 6.14 
at Misoh = 1.2 x 10 14 /i _1 M Q , but as we argue above, 
their mass is underestimated, so this should really cor- 
respond to Migob = 1.5 x 10 14 /i _1 M Q . Rescaling, with 
13 = -0.1, we find their value at M 20 Qb = 10 14 /i _1 M Q is 
C2006 — 6.4 ± 0.3. The central value is in good agreement 
with the value we find when following their procedure 
of correcting the halo center misidentification with mock 
catalogs, c 2 oob = 5.8 ± 0.7. The statistical error is larger 



in our case because we do not use the small scale infor- 
mation. Note that because the mocks are likely not to 
be a completely realistic description of the effect, [45[ at- 
tach a relatively large systematic error of 30% on top of 
the relatively small statistical error they obtain. Instead, 
we trade statistical power for reduction of systematic er- 
ror by using the fits from 0.5/i _1 Mpc, where the effects 
of halo center misidentification are less severe, in which 
case we find C200& = 4.6 ± 0.7 for our standard fit with 
(3 = 0.1 and no offsets, and very similar values also for the 
fits from scales above l/i _1 Mpc, either with or without 
offsets. We conclude that our result is relatively insen- 
sitive to the offsetting procedure, and while the results 
agree within the errors, our concentrations are signifi- 
cantly lower. 

Our results can be compared to those from other meth- 
ods that are used to determine the density profiles. Re- 
cent cluster X-ray and strong-weak lensing studies have 
found that the profile is consistent with the NFW model, 
but in many cases with a concentration that is higher 
than predicted by the concordance cosmology implied 
by WMAP @, I2I [H HH. While some previous work 
concluded that this result implies a higher normalization 
cosmology, this interpretation may be premature. There 
are many alternative explanations that need to be ex- 
plored, such as the use of information from scales below 
100/i" 1 kpc in clusters, where baryons have a significant 
contribution to the density profile and tend to steepen the 
profile, therefore increasing the best-fit concentration. In 
addition, there are significant effects of triaxiality on the 
formation of arcs in the strong lensing regime. In the case 
of X-ray analysis, additional sources of pressure support 
may complicate the reconstruction based on hydrostatic 
equilibrium, since only the thermal pressure can be mea- 
sured directly, while other sources of pressure, such as 
turbulence, cosmic rays or magnetic fields, cannot. In 
some cases, the gas may not be relaxed at all, and the 
hydrostatic equilibrium assumption is invalid. Another 
possibility are selection effects, since the dispersion of 
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the concentration is large and correlates with X-rays or 
strong lensing selection such that only high concentra- 
tion clusters near the mass threshold are in the sample 
26]. It is important to compare the different tracers for 
consistency, and our analysis provides a complementary 
approach that can be compared against these more tra- 
ditional analyses only once these other effects are under- 
stood in more detail. 



D. Comparison to simulations 

Next, we compare our concentration parameter fits to 
theoretical expectations. We show theoretical predic- 
tions on top of the data in Figure El u sing the results 
from the Millenium simulations in [3^| and [93]. Us- 
ing a full sample, including both unrelaxed and relaxed 
halos, [lifind c 20 q £ _ = 4.67[M 20 o c /(10 14 /i- 1 M q )]- 011 
(see also [94j and [95]). First, we move these z = 
results to our mean redshift of 0.22, lowering the am- 
plitude by 1/1.22. Then, carefully converting both the 
mass and the concentration to account for the different 
halo definitions, we find the corresponding relation to be 
c 2 oofc = 7.1[M 2 oofc/(10 14 /i _1 M Q )]^ al . This prediction is 
for the Millenium simulation cosmology with fl™ = 0.25 
and er 8 = 0.9; if we convert to erg = 0.82 [93] by as- 
suming that c oc M^j- 1 0, then the predicted amplitude 
of this relation is reduced from 7.1 to 6.7. Finally, the 
results of [93|, HI] suggest that the concentrations in [39[ 
are too high by 10 per cent at the high masses where we 
have the most statistical power, so this would bring the 
predicted value to 6.0 for the WMAP5 cosmology. Note 
that the scaling with mass no longer holds at the high 
mass end, where concentration becomes constant and is 
given roughly by C2oofc ^5 — 6. Simulations predict that 
this occurs at masses comparable to or slightly higher 
than our highest mass bin, so we will continue to use 
the power-law scaling with mass in our analysis. To be 
more quantitative we plot the predictions for WMAP1 
and WMAP3 cosmologies at z = 0.22 for the A^oob mass 
definition on Fig. O together with the observational con- 
straints. As shown, the results for the lower normaliza- 
tion cosmology are ~ 2er above our measured concentra- 
tion at 10 u /i _1 M Q (fits 2 or 8 in table|TT|. 

Typically the predicted profiles are derived from N- 
body simulations by fitting an NFW or Einasto pro- 
file to the 3d density distribution and averaging the ob- 
tained profile over all halos of a certain mass. There are 
many reasons why our observational procedure may dif- 
fer from this. One is that we have both scatter in the 
mass-concentration relation and deviations from spheric- 
ity, both of which can change the mean profile of 2-d 
AS when compared to the average density in 3-d. We 
also have scatter in the mass-richness relation, so that 
the assumption of a narrow mass distribution may be 
violated. When using the signal for central galaxies in 
the brightest luminosity bin in the simulations from [4t| , 
which incorporate both scatter in the mass-luminosity re- 



lationship and in the concentration-mass relationship, we 
found that the best-fit concentration can be 20 per cent 
lower than if there is no scatter. However, the compar- 
ison to the expected concentration at the corresponding 
mass suggests that even without scatter, the concentra- 
tion fits can be biased by up to 20%, possibly due to 
deviations of the average profile from NFW or Einasto, 
which show up differently in AE(i?) than in the spherical 
radial profile p(r). Another comparison in [58| also found 
that the concentration derived from AE(i?) can be either 
an underestimate (at low masses) relative to what is de- 
rived from the radial profile or an overestimate (at higher 
masses) , but it is not clear how this result should be ap- 
plied to our analysis since that analysis did not attempt 
to mimic our fitting procedure in detail. 

Another uncertainty arises from the predicted values 
for concentration in existing simulations. 99] find values 
that are about 10-20% higher than the values used above. 
It is not clear how worrisome this is, given that it was 
derived from only a handful of preselected massive clus- 
ters. A more concerning issue is the difference between 
relaxed halos versus all halos when fitting for concentra- 
tions. When [39| use relaxed clusters only when fitting 
for concentration, they find about 10% higher values than 
for the full sample. Similarly, [94[ find a large difference 
between the two, with relaxed clusters having typically 
20% higher concentrations at these masses. On the other 
hand, analysis of relaxed clusters in [93| can be compared 
to the full analysis presented in [95| and the latter gives 
only a 5% reduction in the values of concentrations at 
the halo masses around 10 14 ft. _1 M Q . Our sample con- 
sists of most or all halos above a certain mass threshold 
in a given volume, hence we should compare it to the full 
sample, which could bring the observed and predicted 
values into a better agreement. 

The above discussion suggests that there is some the- 
oretical uncertainty in the predicted values of concentra- 
tions, at a level of 20%, that prevents us from making a 
more quantitative comparison to our results. This uncer- 
tainty could be reduced if exactly the same analysis used 
here on the real data is repeated on a large sample of 
simulated clusters, but doing so requires a large library 
of simulated galaxies and clusters from cosmological sim- 
ulations with a volume comparable or larger to that used 
in the actual data analysis (of the order of Gpc 3 ) yet 
with high mass resolution, for which the next generation 
of simulations will be required, and is thus beyond the 
scope of this paper. 



E. Implications for shear-shear lensing 

Our results also have implications for the theoretical 
interpretation of shear-shear lensing. The weak lensing 
power spectrum quantifies galaxy distortions produced 
by lensing, which is sensitive to all matter in the uni- 
verse, and as such it has long been advertised as be- 
ing insensitive to astrophysical uncertainties present in 
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other tracers such as galaxies. However, this argument 
relies on the assumption that the baryonic effects on the 
distribution of total matter can be understo od. Whil e 
earlier estimates found the effect to be small [lOOl llOlj , 
recent work based on simulations with gas cooling finds 
that this assumption may be invalid jl02l Il03l |. They 
find that the baryons cause a significant redistribution of 
matter within a halo, such that cooling and compression 
of baryons towards the center also makes the dark mat- 
ter more concentrated. This redistribution can increase 
the concentration relat ive to the theoretical model pre- 
dictions by up to 40% jl03l Il04| , and suggests that the 
matter profile is significantly redistributed well outside 
the inner region of the halo where gas has been trans- 
formed into stars. 

Observationally, we find no evidence for such an in- 
crease in concentration. We find that concentrations 
are at the lower end of the range predicted by simula- 
tions even for the low normalization cosmology implied 
by WMAP3, and more so for the latest determinations 
of normalization and matter density [96l \9l\ , which are 
already above the observations by 2a. An increase in 
predicted concentrations by 40% would make the discrep- 
ancy more than 5cr, which is ruled out by our analysis. 
If halo center misidentification were considerably worse 
than estimated, then such a concentration enhancement 
would still be possible, but this would then be inconsis- 
tent with our analysis of halo concentration determina- 
tion as a function of inner radius (see table [TXJ) , where 
we find no evidence of a systematic change in concentra- 
tion with radius. It would also be inconsistent with our 
comparison between the analyses with and without dis- 
placements, for which we find the difference to be within 
the statistical error if the analysis excludes information 
below 0.5/i _1 Mpc. Thus, on scales larger than this, ef- 
fects due to baryon cooling are likely to be small, so that 
their effect on the existing shear-shear measurements can 
be neglected. 

In the future, the statistical error of shear-shear auto- 
correlation measurements will be significantly reduced, so 
the baryonic effects may therefore become significant, but 
at the same time, the measurements of the halo profiles 
will improve as well. Thus, the effect can to some extent 
be corrected for by including the differences between the 
observed and predicted galaxy-galaxy and cluster-galaxy 
lensing profiles in the analysis. One approach to do so is 
through t he ha l o model analys is of the dark matter power 
spectrum [Tol, QUI [To3, Il08l |. which has been shown to 
give a good agreement with the simulations [1031 ] . A nec- 
essary requirement, however, is that our understanding of 
BCG halo center displacements improves, either through 
observations or via simulations. 



IV. CONCLUSIONS 

We used a large sample of 170 640 isolated spectro- 
scopic galaxies, 38 236 groups traced by spectroscopic 



LRGs and 13 823 MaxBCG clusters from SDSS, and ap- 
plied a weak lensing analysis to determine their average 
masses and matter density profiles as a function of halo 
mass. We fit the lensing signal to an NFW or Einasto 
profile, excluding small scales to reduce the effects of 
baryons and misidentification of halo centers. The largest 
scale we use in the fits is 3ft- _1 Mpc for clusters, where the 
large scale structure contribution is still small, but we ac- 
count for it using a halo-halo term with a bias predicted 
by simulations. For galaxies, the largest scale we use is 
500/i _1 kpc, and for groups/LRGs, l/i _1 Mpc. 

Fitting the lensing signal to an NFW profile, we find 
that the cluster concentration weakly decreases with 
mass, c = co{M/M )-P, with (3 = 0.13 ± 0.07 in good 
agreement with predictions from simulations. The mean 
concentration at M2006 = \Q 1A h~ 1 M ( r ) is C200& = 4.6 ±0.7 
(z = 0.22). This value should be compared to the pre- 
dicted value ~ 6 for the best fit cosmological models 
(96l [97] | . The measured concentrations are below the pre- 
dictions, although within 2a. We find very little differ- 
ence between NFW and Einasto profile in terms of the 
measured concentration. 

While there appears to be a mild discrepancy between 
the predictions and observational constraints, there are 
significant uncertainties in the theoretical predictions 
that prevent us from robustly concluding whether there 
is a problem. One task for the future is to repeat exactly 
the same analysis as done here on on a representative 
sample of halos from cosmological simulations covering 
the mass range of observed halos. This should be possi- 
ble in the near future as a new generation of large volume 
and high mass resolution N-body simulations becomes 
available, thus allowing for a more accurate calibration 
of the concentration-mass relation than is possible at the 
moment. 

However, to reduce the systematic uncertainty further 
we also need a better understanding of the displacement 
of BCGs from the halo center, which can be achieved 
either through observational studies [56| or through im- 
proved modeling in cosmological simulations [45| . In this 
paper we attempt to minimize it by using information 
outside the central region where we expect the effects 
from baryons in the central galaxy and from misesti- 
mation of the cluster center to be small. We see no 
evidence of systematic contamination in the sense that 
we find consistent results with and without accounting 
for the halo center misestimation, but we cannot com- 
pletely exclude the possibility that there are residual ef- 
fects at a level comparable to or below the statistical 
error. Even if we increase the measured concentrations 
by this amount, they do not exceed the predicted values, 
and thus we see no evidence of an enhancement in con- 
centrations due t o ba ryonic cooling predicted by some 
simulations [l02l Il03l | . This result bodes well for exist- 
ing and future shear-shear weak lensing analyses, in that 
the baryonic effects are likely to be small and confined to 
small scales, and that by comparing theoretical and ob- 
served profiles as done here, these effects can be corrected 
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for. 
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